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Abstract: The in-medium feature of nuclear force which includes both nucleon-nucleon (NN) and hyperon-nucleon 
(AN) interactions impacts the description of single-A hypernuclei. With the alternated mass number or isospin 
of hypernuclei, such effects could be unveiled by analyzing systematical evolution of the bulk and single-particle 
properties. From a density-dependent meson-nucleon/hyperon coupling perspective, a new AN effective interaction 
in the covariant density functional (CDF) theory, namely DD-LZ1-A1, is obtained by fitting the experimental data 
of A separation energies for several single-A hypernuclei. It is then adopted to study the structure and transition 
properties of single-A hypernuclei in Oxygen isotopes, comparing with several selected CDF Lagrangians. Discrepancy 
is observed explicitly in the isospin evolution of Alp spin-orbit splitting with various effective interactions, ascribed to 
their divergence of the meson-hyperon coupling strengths with increasing density. In particular, the density-dependent 
CDFs introduce an extra contribution to enhance the isospin dependence of the splitting, which is originated from 
the rearrangement terms of A self-energies. In addition, the characteristics of hypernuclear radii are studied along 
the isotopic chain. Owing to the impurity effect of A hyperon, a size shrinkage is observed in the matter radii of 
hypernuclei as compared to their cores of normal nuclei, while its magnitude is elucidated further to correlate with 
the incompressibility of nuclear matter. Besides, there exists a sizable model-dependent trend that A hyperon radii 
evolve with the neutron number, which is decided partly by the in-medium NN interactions as well as the core 
polarization effects. 


PACS: 21.80.+a, 13.75.Ev, 21.30.Fe, 21.60.Jz 


1 Introduction 


The discovery of hyperon, particles containing strange quarks, in 1953 sparked strong interest among experimental 
and theoretical physicists [I]. The ability of hyperons to enter the nucleus and form a system of hypernuclei makes 
them sensitive probes for studying the structure and specific nuclear features. The studies on hyperon behavior 
in the nucleus help us to understand the baryon-baryon interaction in nuclear medium and its effects on nuclear 
properties [2] B]. In addition, hyperons are thought to be produced inside neutron stars [4}{6]. The link between 
hypernucleus and neutron star properties benefits our comprehension of the state of matter in extreme environments, 
as well as strangeness-bearing nuclear force at high densities. In recent decades, a wealth of hypernuclear data 
has been generated through induced reactions of meson and electron beams at various radioactive beam facilities, 
including the Japan Proton Accelerator Research Complex (J-PARC) [7], the Thomas Jefferson National Accelerator 
Facility (JLab) [8], and the Facility for Antiproton and Ion Research (FAIR) [9]. These advanced facilities have 
played a pivotal role in advancing our understanding of strangeness in nuclear physics. Notably, single-A hypernuclei 


have been the most extensively studied, with experimental data covering hypernuclei from 3H to %°8Pb in various 


laboratories [2] [3] [0] MI]. 

When A hyperon enters into a nucleus, various phenomena could be observed. For instance, in (Li, it has been 
found that the size of the ĉLi core is smaller compared to the free space °Li nucleus, as suggested by the measurement of 
the y-ray transition probability from £2(5/2+ — 1/2*) in {Li [I2]. In addition, in }C, it is hinted that the A spin-orbit 
splitting is much smaller than the nucleon’s [13]. Recently, the potential for producing neutron-rich hyperfragments 
at high-intensity heavy-ion accelerator facilities is discussed [14] [5]. The directed flow of hypernuclei (3H and 4H) 
just observed at RHIC for the first time in heavy-ion collisions, providing insights into hyperon-nucleon interactions 
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under finite pressure [16]. These advances highlight the promising prospects for investigating hypernuclear structures 
using the forthcoming high-intensity heavy-ion accelerator facility HIAF [8]. To provide accurate predictions for 
these experiments, researchers have performed detailed theoretical work on observables such as hypernuclear binding 
energy [20], spin-orbit splitting [21}{23], hyperon and hypernuclear matter radius [24428]. Overall, these efforts 
aim to provide valuable insights into the behavior of hypernuclei, and to deepen our understanding of the in-medium 
baryon interactions. 

Due to their ability to provide a self-consistent and unified description of almost all nuclei on the nuclear chart, 
both non-relativistic and relativistic mean-field theories are widely used in the calculation of finite nuclei and nu- 
clear matter, and have been extended to describe hypernuclear systems with strange degrees of freedom during the 
development of theoretical models [29146]. As a key model utilized in this work, the relativistic mean-field 
theory has been extensively developed to study hypernulcear properties such as hyperon separation energy [47], 
spin-orbit splitting [23] [39] [48], hyperon halo [49], hypernuclear deformation [25] [50452], cluster structure [53] and 
drip lines [54]. While most theoretical models have primarily emphasized nonlinear self-coupling interactions for 
studying hypernuclei, there has been a recent study that explores the effective interactions for single-A hypernuclei 
within the density-dependent relativistic mean-field (DDRMF) model [55]. With three distinct fitting approaches, 
they propose six new sets of effective AN interactions and uncover a significant linear correlation between the ratios 
R, and R,,, representing scalar and vector coupling strengths, respectively, between these effective AN and NN 
interactions. 

Recently, a new type of density-dependent relativistic mean-field Lagrangian, DD-LZ1, has been proposed, inspired 
by the restoration of pseudo-spin symmetry (PSS) and nuclear medium effects [56]. This new effective Lagrangian 
has produced satisfactory results in describing the properties of nuclear matter and finite nuclei. With unique 
density-dependent form, DD-LZ1 eliminates the spurious shell closures that appeared in previous RMF calculations, 
and reasonably restores the PSS of high orbital angular momentum near the Fermi energy [56]. Applications with 
this new RMF Lagrangian has been performed for several nuclear many-body characteristics, in both finite nuclei 
with mass ranging from light to superheavy, and neutron star properties with density ranging from low to high. 
For instance, a comprehensive macroscopic-microscopic model was developed to evaluate the total energies for even- 
even nuclei with proton numbers ranging from 8 to 110 [57]. Even with the appearance of hyperon [59], larger 
maximum masses of neutron stars could be obtained with DD-LZ1 than with several other RMF parameter sets, 
providing the possibility that the secondary object observed in GW190814 is a neutron star [60H62]. Utilizing the 
Thomas-Fermi approximation, different microscopic structures of nonuniform nuclear matter were calculated for the 
crust of neutron stars and a unified equation of state was established in a vast density range [64]. The different 
density-dependent behaviors of meson-nucleon couplings impact the microscopic structures of neutron star matter 
with DD-LZ1, affect correspondingly the description on various physical processes and evolutions of neutron stars. 

Apart from dealing with the different nuclear medium effects caused by the interactions themselves, the evolution 
of isospin also leads to significant changes in the in-medium effects of hypernuclei, thereby affecting the description of 
their structural properties. In recent years, a series of refined theoretical studies have been conducted on hypernuclei in 
different isotopic chains using various interaction models. For instance, the no-core shell model has been employed to 
investigate the systematic evolution of the ground and excited state energies in the Helium and Lithium hyperisotopes 
[20]. The antisymmetrized molecular dynamics method has been applied to explore the e low-lying level structure of 
hypernuclei in the Beryllium hyperisotopes [65]. The multidimensionally constrained RMF model has been used to 
study the shape evolution of hypernuclei in the Argon hyperisotopes [5I]. The beyond mean-field approach has been 
utilized to discuss the evolution of p-state energies and composition in the Carbon hyperisotopes [23], as well as the 
hyperon halo structures in the Boron and Carbon hyperisotopes [26}|28]. The studies exhibit the significance of isospin 
role in the description of hypernuclear structure. In fact, with the development of hypernuclear spectroscopy, new 
experiments related to hypernuclei have been initiated, such as the planned measurements in the J-PARC project, 
aiming to study the A hyperon binding energies in neutron-rich hyperisotopes of \7*~'*°Sn [66}(67]. These experiments 
will provide crucial information about the properties of hypernuclei associated with various isospin circumstance. 

In view of the essential role of nuclear in-medium effects on hypernuclear structure and their relevance to the 
isotopic evolution, we aim to further expand the density-dependent RMF model to investigate the structure of 
single-A hypernuclei in Oxygen hyperisotopes. First, we will introduce the theoretical framework of the hypernuclear 
RMF approach in Sec. Then, the induced A-nucleon (AN) effective interactions will be determined by fitting A 
separation energies to the experimental data for DD-LZ1 Lagrangian. To give the results and discussion, the influence 
of nuclear in-medium effects will be studied in Sec. on the isospin dependence of hypernuclear bulk properties, 
hyperon spin-orbit splitting and matter/hyperon radius. Finally, a summary will be given in Sec. 


2 DDRMEF approach for spherical single-A hypernuclei 


To describe single-A hypernuclei within the meson-exchanged type of the relativistic mean-field theory, the co- 
variant Lagrangian density serves as the foundation, which is 
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where the index B (B’) represents nucleon N or hyperon A, with its sum 5>, over nucleon N and hyperon A. The 


masses of the baryon and mesons are given by Mg and my (¢ =0,w",p"), while 04”, Re” and F™ are the field 
tensors of vector mesons w”, 0” and photon A”, respectively. The interaction between nucleon (hyperon) and mesons 
(photon) is involved by the Lagrangian %, 
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Here the A hyperon (namely Yp taken as pa), which is charge neutral with isospin zero, only takes part in interactions 
that are spread by isoscalar mesons. The nuclear in-medium effects are introduced phenomenologically via the 
coupling strengths gg (ggn), which use baryon-density dependent functions in density-dependent RMF (DDRMF) 
approach to define the strengths of different meson-baryon (meson-nucleon) couplings [56] [68]. 

The effective Hamiltonian operator for A hypernuclei can be obtained by performing the general Legendre trans- 
formation on the Lagrange density X in Eq. (1), and it can be written as the sum of the kinetic energy operator T 
and the potential energy operator Ŷ,, 
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here x is four-vector (t,a). Correspondingly, we define interaction vertices Y,,(x) for a various of meson (photon)- 
nucleon (hyperon) coupling channels, which for isoscalar ø and w mesons are represented as 


G,p(@)=4+9>n(2), (6a) 
Gip(z)=+gun(x)y". (6b) 


Notably, both nucleons and the A hyperon can contribute to the isoscalar meson fields. However, for the remaining 
isovector mesons and photon fields, it is expected that their interaction vertices solely connect to nucleons since the 
isoscalar and charge-zero nature of A hyperon, 


Gon (2) =+ Gon (x) Y"F, (7a) 
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As the retardation effects could be neglected in the majority of RMF models, the meson (photon) propagators D, 
(Da) read as 
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The baryons field operator wz in the Hamiltonian can be second quantized in the positive-energy space under 
the no-sea approximation as 


ba(z)=)>_ hek e (9) 


Here, f; represents the Dirac spinor, while c; denote the annihilation operators for state i. Accordingly, the energy 
functional E is determined by evaluating the expectation value of the Hamiltonian with respect to a trial Hartree-Fock 
ground state |®,), 
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Then the binding energy of a A hypernucleus is written by 
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where the kinetic energy functional of baryons is shown by Ekin,g. The contributions of the potential energy functional 
from o and w are denoted by the variables Es, and Eu B. Additionally, Epy and Fem. are used to represent the 
contributions from p and A, respectively. The center-of-mass adjustment to the mean-field is represented by the term 
Esm., While Epair takes into account the contribution from nucleon pairing correlations [69]. 

The role of deformation in single-A hypernuclei has been discussed in various density functional models [23} [25] [70] 
[71], which may generate non-negligible effects on the single-particle energies like in Carbon hyperisotopes [23] [25] [71]. 
To describe single-A hypernuclei, in particularly the Oxygen hyperisotopes discussed hereafter, we just restrict the 
RMF approach to the spherical symmetry. Correspondingly, the Dirac spinor f;(a) of the nucleon or hyperon in Eq. 


(9) has the following form: 
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where the index a consists of the set of quantum numbers (nK) = (njl), and Q,.m is the spherical spinor. Meanwhile, 
the propagators can be expanded in terms of spherical Bessel and spherical harmonic functions as 
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where Q = (V, p), and Ry, contains the modified Bessel functions J and K as 
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In the DDRMF approach, the meson-baryon coupling strengths are adopted as a function of baryon density p», 
which are written by 


JoB (Po) =Gon(0)fon(E) or gen (pr) =Gon(O)e*#?*, (16) 
where E= p /po with po the saturation density of nuclear matter, and 
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The free coupling strength at p, = 0 is represented by g,,(0) in the expression above. To keep the variational 
self-consistency between the energy density functional and single-particle properties, the extra terms in baryon self- 
energies, namely the rearrangement terms, will occur due to the density dependence of the coupling strengths. The 
single-particle (nucleon or hyperon) properties can be determined by solving the Dirac equation, 
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Here the self-energies D2 = Xo, +¥s,g composed by the vector and scalar terms. The scalar self-energy Xs, = XZ p, 
and the time component of the vector one has 


o,e (r) =) E$ a(r) +Er(r), (19) 
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where ¢=w,p for nucleons, and ¢ =w for A hyperon. The self-energies of nucleon or hyperon include scalar one 
Xs g and vector one Xo, g, in which the coupling of isoscalar mesons contributes as follows, 
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B' 

Zs s(r)=+ge(r)) / r°dr'gus(r')pe n (r) Roo(r, r). (20b) 
B’ 


Here, p;,z and p, g represent the scalar and baryon density, respectively [69]. Additionally, the rearrangement term 
“pz appears in DDRMF approach, which contain the summation over all baryons for the isoscalar case of 6 =0,w, 
but only over nucleons for the isovector one. For example, the contribution from o—S coupling is shown as 
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3 Results and Discussion 


In recent years, there has been extensive theoretical research on hypernuclei, particularly focusing on the simplest 
single-A hypernuclei, using RMF and RHF theories. In this section, we aim to extend the effective interaction DD-LZ1 
[56], which has been proven to be successful and promising in determining the properties of nuclear structure in both 
bulk and single-particle aspects, to incorporate A hyperon within the framework of RMF model. To give a comparative 
study and illustrate the role of nuclear in-medium effects, the calculations with DD-LZ1 will be accompanied by 
several existing effective AN interactions within CDF models. These interactions have been significantly expanded 
to incorporate the degrees of freedom of the A hyperon and have yielded many successful findings in the study of 
hypernuclear structure and the properties of dense stars. In detail, density-dependent RMF effective interactions DD- 
LZ1 [56], PKDD [69], DD-ME2, TW99, DDV [72], density-dependent RHF (DDRHF) effective interactions PKO1, 
PKO2, PKO3 [69], and nonlinear RMF (NLRMF) effective interactions NL-SH [I9] and PK1 [73] were selected. 
In these CDF functionals, the w-tensor coupling which has been proved to be essential in reducing A’s spin-orbit 
splitting in hypernuclei is ignored. The Dirac equation is solved in a radial box size of R= 20 fm with a step of 
0.1 fm. For open-shell hypernuclei, we employ the BCS method to account for pairing correlations. As the strength 
of hyperon pairing correlations remains uncertain and may become essential in multi-A hypernuclei, our current work 
solely considers pairing correlations between nn and pp pairs by using the finite-range Gogny force D1S [75], see 
Refs. [76}{79] for details. In addition, the blocking effect should be taken into account for the last valence nucleon or 
hyperon, with a detailed description to the Ref. [69]. 


3.1 Density dependence of AN effective interaction 


For the theoretical study of hypernuclear structure, the AN interaction must be determined first. Since the A 
hyperon is an electrically neutral particle with isospin zero, our focus lies on the coupling strengths between the 
isoscalar-scalar @ meson and the isoscalar-vector w meson with the A hyperon. For convenience, we introduce the 
ratio of the coupling strengths between the meson-hyperon and meson-nucleon, gga/ggn. According to the naive 
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Fig. 1. (Color Online) The calculated A separation energies Ba for the single-A hypernuclei with the RMF effective 
interaction DD-LZ1-A1 in comparison with the experimental data taken from Ref. [3] [10]. 


quark model [80], we fix the ratio of the isoscalar-vector meson coupling strength g.a/gun to 0.666, while the ratio of 
the isoscalar-scalar one g,,/g¢n can be obtained by reproducing the A hyperon separation energy Ba experimental 
data for }°O, 4°Ca, and 3°°Pb [3] [10]. In the fitting process, the hyperon is placed in the 15,2 ground state, and the 
By is defined as follows: 


Ba ($Z) = E(*"Z)—-E(4Z), (22) 


Based on the effective interaction DD-LZ1, we finally obtained a new set of AN interaction, namely DD-LZ1-A1, 
after a fitting process of Levenberg-Marquardt minimization. Then, we calculated the A separation energy Ba as 
well as the single-A energy, with hyperon occupying the ground state 1s;/2 or possible excited states with higher 
angular momentum la. For Ba of DD-LZ1-A1, a remarkable agreement with experimental data is found for most 
of hypernuclei, except for 3°Si with significant deformation and Carbon hyperisotopes with light mass, as shown in 
Fig. Actually, more accuracy description to the light-mass Carbon hyperisotopes could be obtained, by limiting 
the mass region of fitting and taking into account the deformation effects [55]. To investigate the deviation in 
describing the structural properties of single-A hypernuclei using different CDF effective interactions, the coupling 
strength of DD-LZ1-A1 in comparison with other selected CDF functionals are listed in Table |1| One could check 
the root-mean-square deviation A for Ba between theoretical calculation and experimental value, which is defined 
by 


1 N 
A= y Bar — Bey). (23) 
To reveal the systematics, we define A, to be the deviation only for 160, 4°Ca, and 7°°Pb, as well as A, that suitable 
for all hypernuclei. 

From Table {1} it can be seen that different CDF theoretical models have good descriptions for 160, 4°Ca and 
“8Pb, and most parameter sets have good consistency for hypernuclear theoretical calculations and experimental 
data over a large mass range from }7C to *°8Pb. In addition, by comparing three different types of CDF effective 
interactions, we can find that when the ratio of isospin scalar-vector meson coupling strength is fixed to the same 
value, the ratio of isospin scalar-scalar meson coupling strength g,,/gon may satisfy certain linear correlations with 
the ratio of isospin scalar-vector meson coupling strength, which has been systematically explored in some works 
[59] [81]. It should be pointed out that the linear correlation of meson-hyperon coupling strength ratios obtained 
in the RMF framework is obviously not suitable for density-dependent RHF models [69]. 


Table 1. The o-A coupling strengths goa/gon fitted for the DDRMF effective interactions DD-LZ1-A1, PKDD-A1 
69], DD-ME2, TW99 and DDV [72], the DDRHF ones PKO1-A1, PKO2-A1 and PKO3-A1 [69], as well as NURMF 
ones NL-SH [I9] and PK1 [73] by minimizing the root-mean-square deviation A; (in MeV) from the experiment 
values of A separation energies of }°O, 4°Ca and 3°°Pb, where the w-A coupling is fixed to be gua/gun = 0.666. 
A2 represent the root-mean-square deviation between the theoretical calculations and experimental values of A 
separation energies for all hypernuclei shown in Fig. 


DD-LZ1-A1 PKDD-Al DD-ME2 TW99 DDV PKO1-Al PKO2-A1l PKO3-Al NL-SH PK1 


Goa Gon 0.615 0.620 0.620 0.617 0.622 0.596 0.591 0.594 0.621 0.618 
Ay 0.319 0.363 0.245 0.375 0.473 0.265 0.260 0.407 0.916 0.519 
A, 1.810 0.734 0.710 0.684 3.460 0.683 0.527 0.881 1.614 1.184 


In DDRMF approach, the in-medium effects of nuclear force are effectively embedded in the density-dependent 
shape of meason-baryon coupling strength, playing the role in the nuclear structure via the equilibrium of nuclear 
dynamics from various coupling channels. In recent years, analysis based on the equilibrium of nuclear in-medium 
dynamics has been applied to clarify the mechanism of the pseudospin symmetry, the shell evolution, the liquid-gas 
phase transition, and hyperon’s spin-orbit splitting in the CDF models [83]. The delicate in-medium 
balance between nuclear attractive and repulsive interactions may be significantly altered by treating the density 
dependence of coupling strength differently, impacting the description of the properties of nuclear matter and finite 
nuclei with different CDF effective interactions. 
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Fig. 2. (Color Online) Meson-hyperon coupling strengths, namely, the isoscalar goa [panel (a)] and gwa [panel (b)], 
as functions of baryonic density p,(fm~*) for the DDRMF effective interactions DD-LZ1-Al, PKDD-A1, DD-ME2, 
TW99 and DDV, the DDRHF ones PKO1-A1, PKO2-A1 and PKO3-A1, as well as NURMF ones NL-SH and PK1. 


To provide a comprehensive understanding of the in-medium equilibrium in hypernuclei, we present the density 
dependence of coupling strengths for selected CDF effective interactions in Fig. Ba) and Fig. [2{b), corresponding 
to the isoscalar-scalar channel g,, and isoscalar-vector one gwa. There are systematic divergences of the meson- 
hyperon coupling strengths with increasing density among density-dependent RMF, density-dependent RHF, and 
nonlinear RMF effective interactions. Notably, the density dependence of gea and gwa is significantly reduced in the 
DDRHF effective interaction compared to the DDRMF effective interaction. This pronounced reduction in density 
dependence also influences the description of single-particle properties in hypernuclei, such as A hyperon spin-orbit 
splitting [69]. Furthermore, in contrast to density-dependent interactions, the NLRMF effective interaction exhibits 
density-independent characteristics for goa and gua. Consequently, when applying these three types of CDF effective 
interactions to single-A hypernuclei, the systematic deviation could take place in describing the isospin dependence 
of the hypernuclear structure. 


3.2 Bulk properties of single-A hypernuclei in Oxygen hyperisotopes 


To focus on the isospin dependence of single-particle properties, we choose the A hypernuclei and their nucleonic 
counterpart in Oxygen (hyper)isotopes as examples, since they usually take the spherical symmetry. To check the 
accuracy of the chosen interactions in describing the properties of finite nuclei, we first calculated the binding energies 
Fg, charge radii R., and matter radii Rm for Oxygen isotopes using the DD-LZ1 effective interaction. We compared 
the theoretical calculations with experimental measurements, which were taken from Refs. [84H86]. From the results 
in Table |2} we can see that the theoretical calculations and experimental measurements are in good agreement for 
both the binding energies Ez and the charge radii R., for the interaction DD-LZ1. It is worth noting that the 
total matter radius R,, of finite nuclei, unlike the charge radius, still has significant uncertainties based on heavy 
ion reaction experiments. The theoretical calculations of R,, reconcile with the experimental measurements with the 
existence of error bars. 


Table 2. Binding energy Eg, charge radii R. and matter radii Rm of normal nuclei ZTN O, calculated by DDRMF 
effective interaction DD-LZ1, compared to the experimental data 86H88]. 


Nucleus Es (MeV) EZP (MeV) R.(fm) REP: (fm) R,, (fm) REP (fm) 
“O -99.699 -98.732 2.766 2.543 
160 -128.215 -127.619 2.752 2.699 2.619 2.57(2) 
180 -140.017 -139.808 2.749 2.773 2.761 2.64(8) 
20 -150.687 -151.371 2.746 2.868 2.71(3) 
20 -160.364 -162.028 2.746 2.955 2.90(5) 
240 -168.802 -168.960 2.761 3.054 3.18(12) 


Furthermore, we summarize in Table |3| the systematics of the occupied energy level of A hyperon, the single- 
particle energies of A hyperon, the total binding energies, the charge radii, and the matter radii of hypernuclei 
in Oxygen hyperisotopes. In order to give possible reference to hypernuclear experiments, we also calculated the 
strength of electric dipole transition B(E1) between the Alp and Als occupation states. The transition strength is 
expressed as 
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Where e, represents the effective charge of the A hyperon. The integration (f|r|i) can be computed using the radial 
wave functions of the initial and final single-A state, see Ref. [26] for details. 

In the framework of relativistic models, Dirac spinors with both upper and lower components could contribute to 
determining the value of B( £1). However, it is checked that the contribution from the lower component is negligible, 
especially for non-charge exchange channel. Therefore, only the contribution from the upper component is preserved 
in current calculations as a simplification. The inclusion of A hyperon causes the so-called impurity effect inside 
hypernuclei [2]. When the A hyperon is filled in the 1s,/2 state, we can see from the comparison of the total matter 
radii in Table [3] and Table [2] that the introduction of hyperon causes a shrinkage effect on the hypernuclei, which 
is approximately 0.06 — 0.13 fm. Compared with the ground-state results, we observe a significant enhancement in 
A root-mean-square radii when hyperon is filled in higher-lying 1p state. This change in the density distribution of 
hyperon due to different level occupations leads to an overall expansion of the hypernuclear matter radii, different 
from the Als case. Additionally, with the increase of neutron filling, both the hyperon radii, matter radii and B( EF) 
show significant isospin dependence, which can be qualitatively explained by the density-dependence of the coupling 
strength. As indicated in Table |3} when A hyperon occupies the 1p state, its density distribution spreads more 
outward than the nucleonic core. As isospin evolves, more neutrons are filled and their attraction to the hyperon 
increases, correspondingly leading to a significant reduction in the hyperon radius. For B(E£,), its value is determined 
not only by the overlap between initial and final states which are sensitive to the neutron number, but also by the 
effective charge. As a result, the B(E,) values enlarge a little from 15O to {O and go down gradually as isospin 
evolves after N =8. 


Table 3. Properties of single-A states in hypernucleus , +40 calculated with the DDRMF effective interaction 
DD-LZI1-A1, including single-particle energies £s.p., binding energies Æg, charge radii Re, hyperon radii Ra, hy- 
pernuclear matter radii Rm and B(E1) value of the transition from the excited (Alp) state to the ground (A1s) 
state. 


Nucleus A(nlj) Es.p.(MeV) Ep(MeV) R.(fm) R,(fm) R.,,(fm) B(E1)(e?fm”) 


151/2 -14.330 -113.245 2.716 2.115 2.458 

150 Ip1/2 -0.413 -99.590 2.772 5.265 2.813 0.095 
Ips/2 -1.582 -101.002 2.760 4.134 2.674 0.119 
1s1/2 -13.086 -140.507 2.704 2.323 2.555 

170 Ip. -1.059 -128.927 2.756 4.609 2.780 0.109 
Ips/2 -2.278 -130.307 2.746 3.963 2.711 0.121 
151/2 -14.170 -153.506 2.699 2.310 2.682 

190 1pı/2 -1.720 -141.540 2.751 4.291 2.861 0.090 
Ips/2 -3.036 -143.003 2.740 3.824 2.815 0.097 
151/2 -15.394 -165.477 2.695 2.295 2.773 

210 1pi/2 -2.463 -153.079 2.744 4.062 2.927 0.075 
Ips/2 -3.890 -154.635 2.733 3.699 2.890 0.079 
151/2 -16.804 -176.670 2.688 2.277 2.829 

230 1pı/2 -3.285 -163.703 2.737 3.882 2.977 0.063 
1ps/2 -4.841 -165.374 2.725 3.582 2.943 0.066 
151/2 -17.634 -185.728 2.723 2.256 2.969 

250 1pı;2 -3.925 -172.669 2.757 3.836 3.079 0.052 
1p3/2 -5.522 -174.326 2.748 3.562 3.055 0.055 


3.3 Isospin dependence of A spin-orbit splitting 


Motivated by the connection between the density-dependent effective interactions of theoretical models and the 
isospin-dependent properties of nuclear structure, the spin-orbit splitting of A hyperon in hypernuclei, as a promising 
observable in current hypernuclear spectroscopy, will be discussed in this subsection with newly developed DD-LZ1- 
A1 and other selected CDF functionals. The A’s spin-orbit splitting is defined by the difference of A single-particle 
energies between a couple of spin partner states, which is 


Wo a 
AF eG = Eja Sty -1/2 — Eja sty +1/2 (25) 


As shown in Fig. the analysis is carried out for A spin partner states 1p in Oxygen hyperisotopes, with the A 
hyperon occupying its ground state. 

In Fig. [3{a), it is seen that the isospin dependence of AE%, is clearly distinguished with the chosen CDF 
functionals. The curves from NLRMF models tend to be stable with increasing neutron number, while for density- 
dependent RMF or RHF functionals the splitting enlarges generally with isospin. Among them, DD-LZ1-A1 exhibits 
the most significant isospin dependence. Besides, it is clear that the smaller A spin-orbit splitting is predicted by 
DDRHF compared to RMF, which has been illustrated as a result in single-particle properties since the dynamical 
equilibrium between nuclear attraction and repulsion is dramatically changed with the appearance of Fock terms 


To better understand the evolution of A spin-orbit splitting with isospin, we could decompose AE 4, into various 
parts according to its source of the kinetic or potential energy. The values are obtained by left-multiplying the 
transferred Dirac spinor to the Dirac equation Eq. (18), and separate the integrated contributions from different 
self-energie terms. For instance, AF,., comes from the contribution of the rearrangement term iz to A self-energy 
Xo a, as seen in Eq. (19), due to the density dependence of meson-hyperon couplings. Consequently, the rest 
one from the kinetic energy and the density-independent potential energies could be summed over, which means 
AE into+w = AES — AE jea, as discussed in Fig. [3{(b). 


AEkin+o+0 


10 12 14 16 
Neutron number 

Fig. 3. (Color Online) The spin-orbit splitting of Alp spin-partner states as a function of neutron number N for the 
ground state in nae © hypernuclei [panel (a)], and its contribution A Fiin+c+w from the sum of the kinetic energy, 
the density-independent potential energies of ø and w channels [panel (b)], as well as the rearrangement terms A Frea 
due to density-dependent meson-hyperon couplings [panel (c)]. The results are extracted from the calculations with 
the DDRMF effective interactions DD-LZ1-A1, PKDD-A1, DD-ME2, TW99 and DDV, the DDRHF ones PKO1-A1, 
PKO2-A1 and PKO3-A1, as well as the NLRMF ones NL-SH and PK1. 


It is observed that the values of A spin-orbit splitting are primarily determined by AEkin+o+4w. However, the 
isospin dependence of the splitting is weakly controlled by AFyiny.4. except for 13O. Attributed to the occupation 
of vlpı;2 orbit, the A spin-orbit splitting predicted by various CDF functionals systematically reduces from }°O 
to 170. As has been illustrated in Ref. [69], the spin-orbit coupling potential of hyperon is determined mainly by 
the radial derivative of the self-energy X^. In general, the more neutrons are filled into hypernuclei, the larger the 
density circumstance where the A hyperon is housing. Thus, if the model is density dependent like DDRMFs and 
DDRHFs given in Fig. the meson-hyperon coupling strength then weakens and AE%, should become smaller 
correspondingly as the neutron number increases. As seen in Fig. [3{b), such a reduction in AFEkin+o+w is remarkable 
from XO to 4’O, and relatively less significant at larger neutron numbers. 

Different from the NLRMF case, the density-dependent CDFs introduce extra contribution to reinforce the 
isospin dependence of the splitting, as demonstrated in Fig. [3{c), which cancels the reduction trend in AFyinyo4u 
overwhelmingly and finally leads to the enhancement of AE{, with increasing neutron number in Fig. [3{a). In fact, 
the contribution AF,., to A spin-orbit splitting is originated from the rearrangement terms of A self-energies Yo, 
which according to Eq. depends on the density slope of the meson-hyperon coupling strength. As the neutron 
number increases, the density scenario where A lives could get more intense, consequently weaker density dependence 
of the meson-hyperon coupling strength, smaller density slope as well as the suppressed value of AF,,,. Therefore, 
the link between the isospin evolution of A spin-orbit splitting and the in-medium behavior of AN interaction with 
baryon density is elucidated from the discussion on Oxygen hyperisotopes. In consequence, possible experimental 
constraints on AE, along the hyperisotopes could assist us further in understanding the in-medium effects of nuclear 
force. 


3.4 Isospin dependence of matter and hyperon radii 


In the properties of hypernuclear structure, not only the A spin-orbit splitting but also the A impurity effect could 
exhibit the information of in-medium nuclear interactions. In Fig. [4{a), we selected DDRMF functionals DD-LZ1-A1 
and DD-ME2, DDRHF’s PKO1-Al1 and NLRMF’s PK1, to illustrate its influence on the matter radii of Oxygen 
(hyper)isotopes, where the solid and dash-dotted lines correspond to the calculated results for single-A hypernuclei 
and their nucleonic counterpart in Oxygen (hyper)isotopes, respectively. The matter radius R,,, in hypernuclei goes 
up monotonically as the neutron number increases, regardless of the specific model used, where a steep leap from 
230 to 30O corresponds to the effect of new occupation in 1281/2. 
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Fig. 4. (Color Online) The variation of the matter radii of hypernuclei [Panel (a)] and their 1°O core [Panel (b)] 
in Oxygen (hyper)isotopes with respect to the neutron number, with the A hyperon occupying the 1s1/2 ground 
state. The solid and dash-dotted lines represent the calculated results for hypernuclei and normal isotopes without 
hyperon, respectively. The results were obtained with the CDF functionals DD-LZ1-Al, DD-ME2, PKO1-A1 and 
PK1. 


Although divergent values given for Oxygen isotopes without hyperon, all of the selected models are getting closer 
in size of matter radii for hypernuclei, implying R,,, of hypernuclei as a possible model-independent observable. It is 
evident that the matter radii of Oxygen hyperisotopes contract as compared to their nucleonic counterparts, namely 
the size shrinkage due to the impurity effect of the A hyperon. However, the shrinkage magnitude appears to be 
strongly model dependent. Among them, the DDRMF effective Lagrangian DD-LZ1-A1 yields the largest difference 
between the solid and dash-dotted lines, whereas the NLRMF one PK1 shows the smallest disparity. By checking 
the bulk properties of nuclear matter within these CDFs, it is verified that the shrinkage magnitude correlates well 
with the incompressibility, which is 230.7 MeV for DD-LZ1, 250.8 MeV for DD-ME2, 250.2 MeV for PKO1, 282.7 
MeV for PK1, respectively [90]. In fact, the larger the incompressibility K is, the harder the nucleus is 
contracted by the exerted attraction from the filled hyperon inside, consequently the weaker size shrinkage effect in 
the calculated matter radii. The similar relation could be found from the Table II of a work on the isoscalar giant 
monopole resonance of hypernuclei, where the effective nuclear incompressibility modulus was extracted [91]. 

To further distinguish the effects of different interactions on the description of hypernuclear structure, we inves- 
tigate the isospin evolution of the A hyperon radius Ra in Oxygen hyperisotopes using all selected CDF effective 
interactions, as shown in Fig. It is seen tangibly that Ra evolve diversely along Oxygen hyperisotopes with 
different CDF effective interactions. Some effective interactions, like PKO3-A1l, DD-ME2, DDV, and DD-LZ1-A1, 
exhibit a reduced Ra with increasing neutron number. Especially, DD-LZ1-A1 gives the smallest hyperon radii a- 
mong all chosen CDFs, and an strong declining trend. In fact, the core polarization effect due to A hyperon plays a 
significant role in this evolution. When A occupies the 1s,,/2 state, its density distribution is concentrated inside the 
hypernucleus. As a result, the A’s coupling or attraction with the nucleons in the core (here corresponding to '°O) 
appears relatively stronger than that with the valence nucleons. Hence, the evolution of the hyperon radius could be 
comprehended more or less by the size change of the core with respect to the neutron number. 

The variation of the matter radii for the '°O core in Oxygen (hyper)isotopes is plotted in Fig. K) with respect 
to the neutron number. From N =8 to 14, in contrast to the situation of total matter radii Rm, there is no consistent 
isospin dependence for the selected CDFs in the core radius R$%™° with increasing neutron number. The nonlinear 
RMF functional PK1 exhibits a significant increasing trend with isospin, while the density-dependent RMF one 
DD-LZ1-A1 shows a noticeable decrease. Consequently, the hyperon radius Ra exhibit a similar isospin dependence 
resulting from the core polarization effect, determined mainly by the various isospin properties of CDF functionals in 
nucleon-nucleon channels. From such analysis, the importance of nuclear in-medium effects in affecting the hyperon 
radii is unveiled. So the divergent isospin evolution of Ra given by the CDFs with different density dependent 
meson-baryon couplings makes it a valuable tool to elucidate the in-medium behavior of nuclear force. 


4 Summary 


In summary, considering the significance of nuclear in-medium effects in nuclear many-body problems, such as 
eliminating the spurious shell closures, we expanded the newly developed DDRMF Lagrangian DD-LZ1 to incorporate 
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Fig. 5. (Color Online) The variation of the A hyperon radius with respect to neutron number in Oxygen hyperisotopes 
is shown. The A hyperon is filled in the 1s1/2 ground state. These results were obtained from calculations using 
all selected CDF effective interactions. 


the A hyperon degree of freedom and determined the AN effective interaction by fitting the experimental data of A 
separation energies for several single-A hypernuclei. Subsequently, with several other CDF functionals, the features 
including A separation energy and B(E1) transition, and the evolution of the spin-orbit splitting as well as the 
characteristic radii were analyzed in detail along the Oxygen (hyper)isotopes. 

By comparing the results obtained from different CDF models, we further investigated the crucial impact of nuclear 
in-medium effects on accurately describing the properties of hyperon, both in terms of their bulk and single-particle 
properties. For the 1p spin-orbit splitting of the A hyperon, significant differences in the isospin dependence are 
observed among the selected CDF effective interactions in Oxygen hyperisotopes. As the neutron number increases, 
the density circumstance where the hyperon is housing gradually increases, which causes the meson-hyperon coupling 
strengths that determine the hypernuclear properties to change as well. In particular, the density-dependent CDF 
effective interactions introduce additional rearrangement terms that significantly enhance the isospin dependence of 
the A spin-orbit splitting, leading to more distinct variation of AEĝ, with neutron number in DDRMF and DDRHF 
models. 

The evolution of the hypernuclear matter radius with isospin was further investigated. Significant model depen- 
dence in the magnitude of size shrinkage due to the inclusion of A hyperon is observed, where the DDRMF functional 
DD-LZ1-A1 displays the largest shrinkage effect. The result was then explained by an anticorrelation between the 
incompressibility coefficients K of nuclear matter and the hyperon radii Ra, providing us a possible way to constrain 
the hyperon distribution inside a hypernucleus from better-determined bulk properties of nuclear matter. Addition- 
ally, it is found that the isospin evolution of the hyperon radius is primarily influenced by the density-dependent 
behavior of the chosen CDF functional in NN interaction channel via the procedure of the core polarization. Thus, 
the sensitivity in depicting these hyperon-relevant properties in CDF models with a various of different meson-baryon 


couplings holds us great potential to elucidate nuclear in-medium nature in both AN and NN channels. 
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